Hawking temperature and higher order calculations 
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Hawking radiation has recently been explained by using solutions of wave equations across black 
hole horizons in a WKB approximation. Higher order calculations using both usual and non-singular 
coordinates are found to change the solution for zero spin, but this change is not an alteration of 
the Hawking temperature. For spin 1/2, there is no correction to the simplest form of the solution. 



A classical black hole has a horizon beyond which nothing can escape, but quantum theory alters the situation as 
shown by the theoretical discovery of radiation from black holes For a Schwarzschild black hole, the radiation, 
which is thermal, has a temperature 

~ A-KTh ~ SttM' '''^^ 

where Vh gives the location of the horizon in standard coordinates and M is the mass of the black hole. This was 
derived by considering quantum massless particles in a Schwarzschild background geometry. The derivation being 
quite complicated, attempts have been made to understand the process of radiation by other methods. In 0, a 
path integral study was made, and analytic continuation in complex time used to provide further evidence for the 
temperature Tjj- Furthermore, the propagator in the Schwarzschild background was shown Q to have a periodicity 
in the imaginary part of time with period 47rr/j = SttM, again suggesting the same temperature. There is also an 
argument involving a conical singularity on passing to imaginary time, which can only be avoided if the standard 
Hawking temperature is chosen. 

Later, other attempts were made to understand the emission of particles across the horizon as a quantum mechanical 
tunnelling process [1, |3| . A WKB treatment was used to reproduce the standard temperature Th @ . Some ambiguities 
were resolved in [g!] and clarified using non-singular coordinates in ^ . Various related discussions can be found in . 
We shall discuss here corrections to the wave function when higher orders in Planck's constant are considered. 

A massless spin zero particle in the Schwarzschild background is described by the Klein-Gordon equation 

4^5p(.9^"V^5,<i>) = 0. (2) 



One writes 

$ = exp(-l^) (3) 

and obtains, to the leading order in Ti, the equation 

g'^'df.Sd^S = 0. (4) 
If we consider spherically symmetric solutions and use separation of variables to write 

S = Et + C + Soir), (5) 
where C is a constant arising from the integration of ^ — E, the equation for Sq becomes 

- + (1 - '})S',{rf = (6) 
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in the Schwarzschild metric. The formal solution of this equation is 

S,{r)^±E f (7) 

The two solutions correspond to the fact that there are incoming and outgoing solutions. There is a singularity at 
the horizon r = rh, which is taken care of by skirting the pole and changing r — rh to r — — ie. This yields 

So{r) ^ ±E[r + Th ■ iTT + rh f drP{^—)], (8) 

J r -Th 

where P() denotes the principal value. Let us write the incoming and outgoing solutions for r > as 

S,n = Et + C + E[r + rh-iTr + rh f drP{^—)], 

J r -rh 

Sout = Et + C-E[r + rh-in + rh [ drP{^—)], (9) 

J r -rh 

The real part of C is arbitrary, but if the imaginary part is ignored, strange consequences follow. The imaginary 
part of Sout implies a decay factor exp{—TrrhE/h) in the amplitude, and a factor exp(—2'KrhE/h) in the probability, 
which is not consistent with the known Hawking temperature. Further, the imaginary part of Sin also implies a decay 
factor exp(— Trr/ji?/?!) in the amplitude, and a factor exp(— 27rr/j£'/?i) in the probability: the motion is in the inward 
direction, so this factor in the probability is not exp{+2iTrhEjh) (which is greater than unity) as might appear at 
first sight. Thus absorption is inhibited, as in a white hole To sort this out, one should use coordinates which 
are smooth across the horizon [7|, but a simple method is to make the imaginary part of C cancel the imaginary 
part of Sin, so as to ensure that the incoming probability is unity in the classical limit - when there is no reflection 
and everything is absorbed - instead of zero or infinity The need for the boundary condition of unit incoming 
probability arises because the Schwarzschild coordinates do not distinguish between a black hole and a white hole, and 
a specific choice has to be made explicitly J"]. Then 

C = -inrhE+iRe C), 

Sout = Et - E[r + rh ■ 2t7r + rh f drP{^—)] + {ReC), (10) 

J r -rh 

implying a decay factor exp(— 27r7'/ii?/?i) in the amplitude, and a factor exjp{—4:T:rhE /h) in the probability, in confor- 
mity with the standard value of the Hawking temperature. 

In the above calculation, S has been treated only to the lowest order in Planck's constant. We shall examine what 
happens in higher orders. For the Schwarzschild metric, the full equation for a spherically symmetric S reads 



[(1 _ ^)-i^2 _ _ '±^g,2^ ^ -ih[r^{l - ^)-'S - (r^il - ^)5')']. (11) 



It may be pointed out that [§] does not include the r^ factors in this equation which come from ^/—g in three space 
dimensions and thus differs from us. We write 

S^Et + C + Soir) + Si{r) + S-iir) + (12) 

where Sn{r) is of n*'' degree in h. The terms of n*'' degree in h for n > yield 

n 

^'(1 - y)T.^rSn-r = '^r'i^ - y)S'n-i]'- (13) 

r=0 

This equation expresses Sn^r)' in terms of lower order S'^^-s and can be solved recursively. In particular, 

S[ = (14) 
r 

S2 = T^. (15) 
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Odd order <S'n(r) are seen to be purely imaginary and of the same sign for outgoing and incoming solutions. Even 
order 5'„(r) are real and differ in sign for the two solutions. 

We may also consider advanced Eddington-Finkelstein coordinates which are non-singular at the horizon. In this 
case, the equation for a spherically symmetric solution reads 

(r^ - rrh)S'^ + 2r^ES' + in[2rE + (r^ - rrh)S" + (2r - r,,)S"] = 0. (16) 

To zeroeth order, one gets two solutions for S", corresponding to incoming and outgoing cases 0: 

^^-0 or (17) 

The incoming solution shows full absorption and the outgoing solution leads to the standard Hawking temperature 
without any need for a C . The next terms are 

These are identical with the corrections found above for Schwarzschild coordinates. In fact, all corrections are as in 
that case. 

As regards the significance of these corrections to S, we first note that Si is not really a quantum correction, because 
Si/h is independent of fi. 

Si/h = ~i\nr, (19) 

meaning that the wave function $ gets a multiplicative factor l/r, which was only to be expected in three space 
dimensions. 5*2, like other even order corrections, is real and does not play any significant role. The first significant 
correction is 

which alters the amplitude of $ for both incoming and outgoing solutions. However, this alteration depends not only 
on E through 1/i?^ but also on r and vanishes at the horizon r = r/i as well as at infinity. It cannot be regarded as 
an alteration of the Hawking temperature. In general, 

Sn « (21) 

which changes the distribution in a complicated fashion. 

It is interesting to see what happens if the Klein Gordon equation is replaced by the Dirac equation. A spin 
connection has to be introduced. In the massless case, two component wave functions are appropriate. There are no 
solutions independent of the angular variables [lo| . but one can take the angular part to be 1/Vsin0. In this simple 
case, the wave function can be written as 

rvsin^ 

the functions F, G satisfy coupled first order equations 

EE = h{l~—)G' 
r 

EG = -nil - —)F', (23) 
r 

which can be solved exactly. The naive incoming solution is 

F = -iG(X'^ie " 1-^, (24) 



* = ^ L_ e-'^'*/^' t,)' : (22) 



while the naive outgoing solution is 



E = iG(xie'' ■' '-^ . (25) 
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Both solutions contain damping factors exp{—TrErh/K) which lead to factors exp{—2TrErh/K) in the probability. To 
make the absorption complete, as in the case of the scalar, a factor of exp(— jC/?i) can be used, whereupon the usual 
emission probability appears. But there are no corrections to the dependence on E or to the Hawking temperature. 
For more involved angular dependence, the equations for F, G have extra terms 1^ and corrections to the above 
solution start from ^2, with Sn cx E^~" once again, changing the distribution in a complicated r-dependent way. But 
one cannot interpret it as a correction to the Hawking temperature. 

To sum up, we have calculated corrections to the spherically symmetric solution of the Klein-Gordon equation 
in Schwarzschild and advanced Eddington-Finkelstein coordinates and found that these corrections, which alter the 
amplitude, do not involve corrections to the Hawking temperature. For the Dirac equation, exact solutions can be 
obtained for a simple angular dependence and in this case there is no correction to the solution, while other solutions 
get altered, but again without any change in the Hawking temperature. 
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